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Elastic scattering of pions from finite nuclei is investigated utilizing a contem- 
porary, momentum-space first-order optical potential combined with microscopic 
estimates of second-order corrections. The calculation of the first-order potential 
includes: (1) full Fermi-averaging integration including both the delta propagation 
and the intrinsic nonlocalities in the tt-N amplitude, (2) fully covariant kinemat- 
ics, (3) use of invariant amplitudes which do not contain kinematic singularities, and 
(4) a finite-range off-shell pion-nucleon model which contains the nucleon pole term. 
The effect of the delta-nucleus interaction is included via the mean spectral-energy 
approximation. It is demonstrated that this produces a convergent perturbation 
theory in which the Pauli corrections (here treated as a second-order term) cancel 
remarkably against the pion true absorption terms. Parameter-free results, including 
the delta-nucleus shell-model potential, Pauli corrections, pion true absorption, and 
short-range correlations are presented. 

Pion-scattering measurements, in combination with phenomenological descriptions Jl 
the propagation of the pion and the delta in the nuclear medium, have proved useful 



for probing details of nuclear structure. The situation is not, however, entirely satisfactory 
because some of the parameters in these phenomenological descriptions have not been derived 
quantitatively, even though their physical origin is believed to be understood. Chief among 
these is a shift || in the energy of the two-body, pion-nucleon scattering amplitude, which 
is evaluated somewhat arbitrarily. Therefore, in this work we want to see how far we can 
go in describing the dynamics of the pion and the delta starting from a purely microscopic 
approach in which the dynamics (including the energy at which the in-medium two-body 
amplitude is to be evaluated) are completely determined from theory. Such an understanding 
is needed before one can envision making a reliable extension of the theory to higher nuclear 
densities and high temperatures, where the propagation of the pion in the nuclear medium 
plays an important role in both heavy-ion reactions and in astrophysical problems. 

From the results of such a microscopic approach we hope to learn the extent to which the 
existing phenomenologies are in quantitative agreement with the dynamics as understood in 
a variety of contexts, including what is known about the reactive content of the interaction 
(true absorption, quasi-elastic scattering, and correlation effects), delta-nucleus dynamics 
(the delta-nucleus interaction, delta propagation, and the Pauli principle), and the inter- 
play of the reaction dynamics with nuclear-structure effects. Although some calculations 
HfH|§,0 of pion scattering do include higher-order terms coming from these effects, a mod- 
ern, microscopic test of pion-nucleus dynamics that makes contact with all this information 
does not yet exist. 

Such a test of pion-nucleus dynamics must deal carefully with several well-appreciated 
but technically awkward aspects of the dynamics. One is Fermi averaging, which is ex- 
pressed as a complicated three-dimensional integral of the off-shell pion-nucleon scattering 
amplitude over the nuclear density matrix. This integration, when performed without any 
approximations, incorporates exactly both the propagation of the delta and the intrinsic 
nonlocalities that are inherent to a two-body resonating amplitude. Another is the Lorentz- 
covariant kinematics. Finally, the pion-nucleon amplitude utilized should contain explicitly 
the nucleon pole. The neglect of this singularity in the two-body amplitude leads || to an 
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artificially low momentum cutoff that produces a geometrical change in the effective radius 
H of the nucleus. We here make a test of pion-nucleus dynamics within the framework of 
the optical potential which incorporates all of these features. 

The isobar-hole model 0, which was a successful semi-microscopic approach to the 
dynamics, has served as a phenomenological tool to fit various pion- (and photon-) induced 
reactions, including the true-absorption and quasi-elastic channels. Much has been learned 
about pion and delta dynamics from this model. Even more has been learned from the 
abundance of high-precision data that have been taken at the meson factories during the 
ten years since the inception of the model. Our work relies on this progress to generate a 
parameter-free microscopic theory, which we will compare here to elastic-scattering data 
from 80 to 226 MeV. 

The improvements that we feel to be needed are naturally incorporated by working in 
momentum space, where the various amplitudes can be written in relatively simple analytic 
form. One technical advance which is particularly suited to momentum space is |I(J the 
use of "relativistic, three-body, recoupling coefficients." These incorporate exactly Lorentz 
covariant kinematics (including Wigner spin precession), and they provide natural 

variables for performing the Fermi-averaging integral. The first-order optical potential is 
given in terms of them by 
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Here, (fc n /c J 4_i|-?/>£ 4 ) is the target wave function (labeled by its eigenvalues a), which is 
proportional to a momentum-conserving delta function and is a function of the relative 
momentum of the struck nucleon and the momentum of the A — 1 remaining core nucleons. 
The pion-nucleon amplitude {k'^k'AtiyV^k^k^} also contains a momentum-conserving delta 
function and is a function of the relative momentum between the pion and nucleon. The three 
implicit momentum-conserving delta functions that are a part of the matrix elements on 
the right of Eq. [l] eliminate all but one three-dimensional integration (the Fermi-averaging 



integration, which must be performed numerically); one overall delta function is left over, 
and this conserves the total momentum . In our work, t(W) is the free-space pion-nucleon 
T-matrix, but the energies that appear in it are shifted by an amount that is calculated by 
a well-defined prescription designed to minimize the effect of the higher-order terms in the 
optical potential. The kinematics involved in Eq. [I] are those of a relativistic three-body 
problem with momenta k n , k n , and k^-i] the details of how relativistic recoupling coefficients 



allow one to calculate Eq. [y can be found in Ref. ||13|| . We work with invariant amplitudes 



TUH that are free of kinematic singularities and utilize invariantly-normed wave functions; 



these introduce phase-space factors into the calculation, which are also treated exactly by 
working in momentum space. 

From our discussion, it is evident that in a momentum-space approach the lowest- 
order optical potential as we formulate it is quite general and can be evaluated without 
approximation. In this sense, our work improves not only the phenomenological optical 
model but also on numerous aspects of the isobar-hole model [[7], which were both 
expressed in coordinate space, where nonlocalities are not as easily handled. The propagation 
of the delta was fully incorporated in the isobar-hole model, but the integration over the 
nonlocalities associated with the two-body amplitude were approximated by factorization- 
an approximation that necessitates a nonnegligible correction, particularly for lighter 



nuclei. To deal with Lorentz-covariant kinematics, expansions and further factorizations of 
integrals were made. Additionally, the pole in the two-body amplitude was neglected, a 
choice which we have been particularly careful to avoid in order to eliminate the possibility 
of a spurious geometrical change in the effective radius of the nucleus. 

Given that we are able to calculate the first-order optical potential without approxi- 
mation, there remains the question of how to organize many-body theory (in particular, 
choosing the energies of the nucleon and the delta in the medium) to optimize its rate of 
convergence. The role of the energy W a in Eq. [L] is quite important in this regard because 
the half-width of the delta resonance (55 MeV) is the same size as typical energies that 
characterize nuclei. Thus, the results of a calculation will be very sensitive to how the ener- 



gies that constitute W a are chosen. W a is defined covariantly as the energy available in the 
center-of-momentum frame of the pion-nucleon system, 

Wl = Wl n -(K + k n ) 2 , (2) 

with W nn defined as the energy available to the tcN pair in the pion-target center-of- 
momentum frame, 



W nn = W - \j(k„ + k n ) 2 + m 2 A _ x , (3) 

and Wq = S, the invariant square energy of the reaction. The mass of the A — 1 system, 
tua-x, differs from the mass of the A-body target, tua, by a nucleon mass and a binding 
energy, = tra-i + rn n + Ef,. In their nonrelativistic limit, Eq. || is known |16| as the 
"three-body energy denominator." 

Utilizing the definition of W a given in Eq. ^| produces needlessly large higher-order 
corrections [fT~6| , p~Tf] in the many-body expansion. This is because the delta-nucleus shell- 
model potential, Ua, which is generally believed to be nearly equal to the potential energy 
of a nucleon in the nucleus, has not yet been included. Including the effects of Ua in the 
T-matrix causes an effective downward shift in the position of the resonance that tends to 
cancel the upward shift caused by the nucleon binding energy in Eq. [j[ To incorporate 
this effect, we have proposed |L7| a treatment of W a in Eq. |^ that includes the Ua in a 



first approximation via an energy-dependent and target-dependent energy shift. This shift, 
called the mean spectral energy, E ms , is derived in Ref. |L7j and may be calculated by 



E (w /rf 3 r0(-)*(r)0W(r)p(r)[/ A (r) 
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where Ua{t) is taken to be equal to the shell-model potential of a nucleon. In Fig. [I] 
we present results for 7r + elastic scattering from 12 C at 80, 100, 148, 162, and 226 MeV. 



The data are from Ref. ||T8[ The dashed line is the result of using the full lowest-order 
optical potential, including E ms . The effects including this shift together with the binding 
of the struck nucleons, is quite substantial jT7| . At all energies shown here we find that the 



inclusion of Ua is not only significant but moves the results remarkably close to the data. 



At this point, the agreement of the theoretical results with the experimental data is 
surprising, because there remains much that has not been considered. We know, for example, 



that the pion true-absorption channel is about one-half [H| of the total reaction cross 
section. The Pauli principle [P,^|,pO[ also should play a significant role in the scattering of 
the light-mass pion from the heavier nucleon. The p-wave character of the pion-nucleon 



interaction produces nonnegligible correlation corrections which enter in the form pi] of the 
Ericson-Ericson-Lorentz-Lorenz correction. We will next include each of these higher-order 
terms. The results will provide a test of our understanding of each piece of the physics and 
the role that it plays in pion-nucleus dynamics. 

In order to utilize existing calculations of the second-order terms, we will make extensive 
use of the local density approximation. For the Pauli and true-absorption terms, we utilize 
the functional form of the second-order corrections as derived in Ref. [p]], 

U {2 \k',k) = X^k-k' p^(k-k f ) , (5) 

where p^ is the Fourier transform of the square of the target density. Microscopic calcula- 

(2) 

tions of higher-order terms yield a coefficient A which itself depends weakly on r. In the 
same spirit as the mean spectral energy calculation, we may define the radius R2 at which 
the pion interacts in a finite nucleus by 

R ^ = Jd 3 r 0^)*(r)(^(r)p(r)r 
/ d 3 r <f)^*(r) T (r) p(r) 

In Table | we give the value of R 2 and the density p(R2) calculated for various pion energies 
for 12 C. We note that over this energy region (80 MeV < T v < 315 MeV) the interaction 
is confined to the nuclear surface and low densities. Here po = 0.16 fm -3 (nuclear matter 
density) and the pion distorted waves are taken from Ref. JI|]. 

The Pauli exchange term can be taken directly from Ref. |3|] evaluated at the density 
p(i?2)- We extend the term by including rho-meson propagation in the intermediate state. 
We omit pion distortions for the intermediate pion to avoid including multiple reflection 
corrections in the Pauli term. The Ag coefficients are given in Table |. The dotted curve 
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in Fig. p] gives differential cross sections resulting from adding the second-order Pauli 
correction to the lowest-order calculation. We see that the Pauli correction is (a) large and 
(b) completely destroys the nearly-quantitative agreement of the dashed curve. 

We will include true absorption by introducing a A parameter determined from the 
spreading potential of the delta-hole model. These two terms cannot be equated directly 
because the the spreading potential occurs in the denominator of the delta propagator. We 
can make the correspondence by first isolating the P33 partial-wave contribution to the 
lowest-order optical potential and expressing it in a resonant form. The difference between 
this potential evaluated twice, once with width Tq + Im W sp and then with width Tq (the 

free width), is a true-absorption potential that can be expanded at low density to give a 

(2) (2) 
Aq independent of r. Rather than expanding, however, we determine Ag by matching this 

difference to Eq. [|at the radius R 2 - The resulting values of are given in Table |. We see 

(2) (2) 

that, at all energies, there is a large cancellation between A (Pauli) and A (spreading), 
yielding a small total second-order correction. The solid curve in Fig. [1] gives the differential 
cross sections obtained when E ms , X^ (Pauli), and Aq 2 '* (spreading) are all included. The 
cancellation of the Pauli and spreading terms is evident. 

Finally, we also include the correlation (or LLEE) corrections. It has been shown W2 
that the LLEE effect can be included in the delta self-energy by a modification 

27 V m n 



oE ms = — p , (7) 



where £ is the usual Lorentz-Lorenz parameter. The value of £ depends on the range of 
the short-range repulsive correlations between nucleons, the range of the pion-nucleon form 
factor, and the strength of the delta-nucleon interaction. We will allow for some uncertainty 
in the LLEE-parameter £. The minimum value that is reasonable is about £/3 = 0.08, which 
results from a pion-nucleon monopole cutoff of 800 MeV/c and no AN interaction. The 
maximum value of £/3 is 0.23, which arises from a cutoff of 990 MeV/c and includes a 
AiV-interaction contribution. This is a value that would give the real part of the delta- 
hole spreading interaction, which corresponds to SE ms = 23 MeV. The final result of this 



work is given by the shaded area between the solid curves in Fig. ^| (corresponding to the 
range 0.08 < £/3 < 23). These results combine the first-order potential in which the delta- 
nucleus potential is included via the mean spectral energy with Pauli, true-absorption, and 
correlation corrections. 

The agreement with the data shown in Fig. Q is not exact, but it is remarkably good for 
a parameter-free calculation. Discrepancies could be due to the fact that our treatment of 
the second-order corrections is neither exact nor totally consistent (we have taken the true- 
absorption term from the delta-hole model). For these reasons, it is probably unwarranted 
to conclude that the smaller value of £/3 = 0.08 is preferred, even though this result is 
everywhere closer to the data. Firm conclusions should await a more thorough, internally 



consistent treatment p3] of all the higher-order terms. We are motivated to pursue this 
treatment because our present calculation is intriguingly close to the data. 

We have for the first time combined a contemporary momentum-space calculation of the 
first -order optical potential with microscopic predictions of the effects of the delta-nucleus 
interaction, Pauli corrections, pion true absorption, and short-range correlations. We have 
seen that convergence of the expansions appears to be enhanced throughout the resonance 
region by (1) collecting (via the mean spectral energy approximation) together with 
binding corrections into the first-order optical potential, and (2) collecting the Pauli and 
true-absorption terms together. This result supports our perturbative approach |23[ to 
calculating the optical potential. Finally, the good results that we find from 80 to 226 MeV 
with no adjustable parameters suggest that pursuing calculations of greater accuracy for the 
second-order terms might yield a definitive determination of the short-range correlations 
(i.e., the parameter £) and the delta-nucleus interaction, Ua- 
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FIGURES 

FIG. 1. The differential cross section for elastic scattering of 7r + from 12 C at the energies 
indicated on the figure. The dashed curves are a complete lowest-order optical-model calcula- 
tion including the delta-nucleus interaction through the mean spectral approximation. Dotted: 
the complete lowest-order optical potential and the second-order Pauli corrections are included; 
solid: the complete lowest-order optical potential, the second-order Pauli, and the second-order 
spreading potential are all included. The data are from Ref. jjD|. 

FIG. 2. The same as Fig. |l|, except the shaded area includes the full lowest-order optical 
potential, Pauli and spreading corrections, and the LLEE correlation corrections. The two curves 
forming the boundary result from the LLEE parameter £ equal to .08 (the lowest curve in the 
forward direction) and equal to 0.23. 
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TABLES 



TABLE I. Parameters for the Pauli and spreading interaction. The A as defined in Eq. || is 
given as a function of the pion kinetic energy T n (Mev) and corresponds to the the density region 
centered about the radius R2 (fm) in 12 C. The units for Aq 2 ^ are fm 3 . 



T^ 


R 2 


P (-Ra)/po 


\ { o ] (Pauli) 


Aq 2 ^ (spread) 


a[, 2) (Sum) 


80 


2.40 


0.289 


-0.40, -1.46 


-0.93, 2.02 


-1.33, 0.56 


100 


2.52 


0.252 


0.08, -1.88 


-1.43, 2.09 


-1.36, 0.21 


148 


2.80 


0.175 


2.50, -1.74 


-3.29, 1.18 


-0.80, -0.56 


162 


2.86 


0.160 


3.20, -0.90 


-3.70, 0.35 


-0.50, -0.55 


230 


2.90 


0.151 


0.54, 2.50 


-1.02, -2.80 


-0.47, -0.31 


315 


2.67 


0.209 


-0.60, 0.26 


0.67, -0.68 


0.07, -0.42 
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